Thermalizing Quantum Machines: Dissipation and Entanglement 
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We study the relaxation of a quantum system towards the thermal equilibrium using tools devel- 
oped within the context of quantum information theory. We consider a model in which the system 
is a qubit, and reaches equilibrium after several successive two-qubit interactions (thermalizing ma- 
chines) with qubits of a reservoir. We characterize completely the family of thermalizing machines. 
The model shows a tight link between dissipation, fluctuations, and the maximal entanglement that 
can be generated by the machines. The interplay of quantum and classical information processes 
that give rise to practical irreversibility is discussed. 



The hypothesis of the quantum appeared suddenly in 
physics as an offspring of thermodynamics, due to the 
work of Planck on the blackbody radiation. In its early 
days however, the new theory developed rather as a form 
of reversible mechanics. One century after Planck's in- 
tuition, the link between quantum mechanics (QM) and 
thermodynamics has been discussed by several scientists, 
and is still an actual field of research . In parallel to 
fundamental issues, the concept of quantum machines has 
arisen recently in the field of quantum information pro- 
cessing Looking back again to history, we sec that 
thermodynamics was born to describe engines. It is thus 
natural to ask whether there is a "thermodynamics" of 
quantum machines, and whether the modern standpoint 
of quantum information can cast some new light on the 
foundations of thermodynamics. After some pioneering 
works j^], these ideas have stimulated many investiga- 
tions in the last months Q. 

In this Letter, we focus on the process of thermaliza- 
tion, that is, the relaxation towards the thermal equilib- 
rium of a system in contact with a huge reservoir (bath). 
More precisely, let ps be the thermal state of the bath 
p a generic state of the system, and p'^ the state of 
the system at thermal equilibrium. A thermalization pro- 
cess is defined by these two requirements: (I) The state 
p'^ ® Pb is stationary; (II) If the system is prepared in a 
state p ^ p^ ^ at the end of the process we have a total 
state psB such that Tyb[psb] — p"" and Trs[pss] - Pb-, 
where Trs,s are the partial trace over the bath and the 
system respectively. 

Thermalizing quantum channels can be realized by let- 
ting the system undergo interactions with the bath that 
are localized in time. Such models, known as collision 
models are admittedly rather artificial as models for 
dissipative processes 0, but are most natural in the 
context of quantum information |8| . The system passes 
through several identical machines U (figure |l]), or sev- 
eral time through the same machine; at each passage, 
it becomes entangled (that is, it shares a part of the 
information encoded in the state) with an ancilla, i.e. 
some degrees of freedom of the bath. At the output of 
the machine, the ancilla is discarded into the bath: the 
information present in the system has undergone some 



degradation, that depends on the state of the bath and 
on the machine. 
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FIG. 1. The quantum channel: a repeated application of a 
unitary U (quantum machine), that couples the state of the 
system with the state of the bath. 



Our main goal is to quantify the role of entanglement 
in this thermalization process. Since a computable mea- 
sure for entanglement of mixed states is known only for 
states of two two-dimensional quantum systems (qubits) 
we consider a thermalization process in which both 
the system and the ancillas are qubits. Before discussing 
entanglement, we give the family of all the thermaliz- 
ing machines U acting on two qubits, and a fluctuation- 
dissipation theorem for the thermalizing channel that 
these U define. 

The model. We start with a description of the model: 
(i) The system is a qubit, and the bath is a reservoir 
composed of an arbitrary large number N of qubits. The 
free hamiltonian for the whole system is 



N 



Ho = Hs + HB = h[S]+Y,h[i] 



(1) 



where h[k] is the operator acting as h ~ ^Ea^ on the 
qubit k and trivially on the other qubits. The bath 
is supposed to be initially in the thermal state ps = 



I3He 
t 



e 

\{ 
Pi = 
thus 



ta,-Qh{l3E)a^), and P ^ Let Pq = |0)(0 
|1)(1| be the projectors on the eigenstates of a 



and 



C = pPo + qPi , q = 1 - p 



(2) 



with p= ^{1 + tanh(/3£:)). We set > 0, so that |0) is 
the ground state, and p = 1 corresponds to T = 0. 
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(ii) The machine C/ is a unitary operation on C (g) C". 
This means that at fixed time the system interacts with 
just a single qubit taken out of the bath. In our model 
we consider that a qubit of the bath undergoes at most 
one interaction with the system. This assumption is jus- 
tified by the fact that the bath is assumed to very large 
(i.e. "infinite"). Therefore the input state of the ancilla 
is always ^, and we write 



pik+i) ^ Trs [U ip^") ® U^] = [p'"'] 



(3) 



Thermalizing machines. For the model just intro- 
duced, the two requirements I and II read |lO|| : 



Req. I: (^ 



Req. II: 



pPo 



qPi ,Vp (4) 
(5) 



The subscript z is meant to remind that we allow the 
machine U to depend on the eigenbasis of ^, that is on h. 
Conversely, we want these requirements to hold for all p, 
that is for all temperature. 

It is important to notice here the existence of an 
equivalence class. Let u{x) = Pq + e*^Pi, and suppose 
that Uz satisfies (^ and defines a channel T^. Then 
C/^ = (l u{a)) Uz (l (8) w(/3)) satisfies (||) as well, and 
defines the same channel = T^. This is easy to see by 
noticing that u{x)£^u{xy = ^ for all x. This equivalence 
is a consequence of the freedom of choosing the global 
phases of |0) and |1) for qubits in the bath. Having no- 
ticed this, we can proceed to find all the thermalizing 
machines. 

Take first Requirement I: condition (^ implies that 
the subspaces Pq <8) Pq , Pi (8) Pi and Pq Pi -I- Pi Po 
must be invariant under the action of U. In fact, on 
the l.h.s. UzPq 'Si PqU^ appears with the weight p^, 
Uz (Po ^Pi+Pi® Po) with the weight p{l —p), and 
Uz Pi®Pi U^ with the weight (1— p)^. Since we want con- 
dition (^) to hold for all p, the three subspaces must be 
separately invariant. This implies [Uz,Ho] = 0: the sum 
of the one-qubit energies is conserved by the interactions. 
By inspection, one can see that, up to a global phase fac- 
tor, the most general unitary operation that leave these 
subspaces invariant is parametrized by five angles; only 
three of them are left if we choose a suitable represen- 
tative element in the equivalence class discussed above. 
Precisely, all unitary operations that fulfill the condition 
(0) can be chosen of the form 



Uz((b,e,a) 



|0)|0) 

|0)|1) 
|1)|0) 



|0)|0) 



(c|0)|l) 
(c|l)|0) 



|1)|0)) 
|0)|1)), 



(6) 



with c = cos (p and s — sin 0, G [0, ^] , and 9,a ^ [0, 2tt] . 
We turn now to demonstrate that almost all these ma- 
chines satisfy the condition (||) as well. To do this, let's 
write the state of the system after n steps as 

pin) ^d^n)p^ _^ (l-d("))Pl + fc(")|0)(l| + h.C. (7) 



Inserting the explicit form (|[) for Uz into (g) , we find that 
the effect of the map is given by = d'^^^c^ +ps^ 

and fc("+i) = cAfc(") with 



A 



(pe 



(8) 



A straightforward iteration gives and fc(") as a func- 
tion of the parameters S-^^ and fc*^"^ of the input state: 

= (1 - (cos 0)2") p + (cos 0)2" d(o) , (9) 
fc(")=A:(")(Acos0)". (10) 

Thus, whenever 0^0, the iteration of T"^ yields c?*^"^ p 
and fc(") ^ since |A| < 1; that is, almost all 

the machines of the form (^ satisfy Requirement II as 
well. In conclusion, the family of thermalizing machines 
for ~ is composed by the Uz{(t>,0,a) given by ^ 
with0y^O|l|. 

Dynamical equivalence of machines. The dynamics (^) 
of the diagonal term d^"-*, that is the dissipation^ is de- 
termined only by 0. The other parameters 9 and a enter 
only the dynamics ([l0| ) of the off-diagonal term k^'^\ the 
decoherence, through the complex number A given in (|^). 
Actually, a plays a trivial role: it simply redefines at each 
iteration the axes x and y in the plane orthogonal to z. 
Apart from this global rotation, the effect of A can be 
visualized as follows: the state of the system undergoes a 
rotation of —9 (resp., +9) in the (x, ?/)-plane if it interacts 
with the state |0) (resp., |1)) of the bath, which happens 
with probability p (resp., q). This dephasing process con- 
tributes only to the decoherence rate. Note also that |A| 
is unchanged if one replaces 9 hy tt + 9. Guided by these 
considerations, we say that two machines Uz{4>, 9, a) and 
Uz{(j),9 + m:,a') that differ only on the value of a are 
dynamically equivalent. We choose 

Vzi^,9)=Uz{q^,9,0) = Uz{(l>,9,a) [u{a)^u{-a)] (11) 

as representative element of the class. The Vz{4>,9) are 
diagonal in the Bell basis: 



Vz 



Poo + Pii + e'(''+'^^P*+ + e^^^-^)P*- (12) 



with = ^(|01) ± |10)). The hamiltonian represen- 
tation is easily derived: 14(0, 9) = e*^ QtHi4>,e} with 



Hi(l>, 0) = -[(!) {a^ (^a^ + ay(E)(Ty)-9az(i) dz] • (13) 

Finally note that we can handle the dissipation and 
the dephasing processes separately, since Vz{(j),9) = 
K(0,O)K(O,0) = VziO,9)Vzicl>,0). By the way, ^,(0,0) 
is a realization of the two-qubit copying machine pro- 
posed by Niu and Griffiths, that defines Eve's optimal 
individual attack on the four-state protocol of quantum 
crypthography [|ll. 



2 



The partial swap. During the whole construction of the 
thermalizing machines, we insisted on the fact that the 
machine may depend on the direction z defined by the 
local hamiltonian h. A natural question is whether any 
of the machines Uz{4>,d,a) is actually independent of z: 
such a machine would thermalize the state of the system 
for all one-qubit hamiltonians h = —E n ■ a. It turns out 
that there exist a unique machine with this property, 
which is 1^(0, —(t>)- This machine is a partial swap, since 



V{(j), -(f)) = e^"^ ( cos 1 + isincj) Usw) 



(14) 



where J/.. 



^(f'~f) ^^"^ swap operation, i.e. it 



is the unitary operation whose action is ® \'4>2) — > 

|V'2)'8)|V'i) for all iV'i) , IV'2) G C . The partial swap con- 
veys the intuitive idea, that at each collision part of the 
information contained in the state of the system is trans- 
ferred into the bath. This machine is the cornerstone of 
the quantum-information process called homogenization 
p^ . For the dissipation, i.e. apart from phase fluctu- 
ations, all thermalizing machines are equivalent to the 
partial swap: V{(j), 6) = 1/(0, 9 + <j))V((j), -0). 

This concludes the characterization of the family of the 
thermalizing machines. In the remaining of the paper, we 
study their properties, first in terms of thermodynamics, 
then from the standpoint of quantum information 

Relaxation times. We'd like to pass from the dis- 
crete dynamics indexed by n to a continuous-time dy- 
namics with parameter t. To perform the limit, we 
set n = t/To, and we let the interaction time tq go to 
zero together with cj) and 9, keeping constant the dis- 



sipation rate 



and the phase fluctuation rate 



TO Ti 

^ = We have (cos^ ,/))" « (1 - 02)^ 



and(|A| COS0)" « [{I -2pqe'^){l - ^)\ 



e ^2 with 



11 1 1 / 9'^\ 

— = — = — l + 4pg lim • (15) 

Thus in the continuous-time limit, the processes of dissi- 
pation (||) and dccohercncc ( pof ) are exactly exponential: 



dit) = e-*/^irf(0) + (1 - e-*/^^i)p. 



\k\it) 



-t/T2 



|fc|(0) 



(16) 
(17) 



with the relaxation times Ti and T2 defined in the usual 
way [pi. For 6* = or at zero temperature, the bound 
Ti > ^2 (see e.g. [§, p. 120) is saturated. 

Fluctuation-dissipation (FD) theorem. A FD theorem 
links the fluctuations at equilibrium and the mechanisms 
of dissipation. Usually, this link is derived in a differ- 
ent framework [Q: in particular, a continuous spectral 
density of the bath is normally assumed, while here Hb 
exhibits finite gaps. Nevertheless, one can define a mea- 
surable quantity associated to fluctuations and dissipa- 
tion. Consider the following protocol. First, the sys- 
tem is prepared in the equilibrium state ^ and is mea- 
sured in the basis of its eigenstates Pq and Pi. Obvi- 
ously, the mean values of one-qubit observables A are 



unaffected by this measurement. Then we let the sys- 
tem undergo n interactions with the bath qubits: from 
the state Pj [j — 0, 1) in which it had been found 
by the measurement, the system evolves into the state 
p(n) ^ xn^p.^^ = (1 - c^")^ c2"P, according to (|). By 

the definition of equilibrium, pPo"''+9Pi"^ = in particu- 
lar, the mean value of A holds unchanged. However, due 
to the information gained through the measurement, now 
we have also access to the following statistical quantity: 



F^-^ = Jp[Tr{5^^^A)Y + q[Tr[5t^A)\ 



(18) 



where 5 



(n) 



(n) 



Pi is the deviation from the mea- 



sured state Pj after n interactions. Fj^ 



(n) 



is a measure of 



the fluctuations of A; the dissipative element can be seen 
through the fact that if F\ ^ 0, then the fluctuations 
have partly erased the information that we had obtained 
through the measurement. Writing Z?^") = (1 — (cosi/))^") 
we find 



1 



2 cosh(/3£;) 



|Tr(^a, 



(19) 



In the continuous time limit, is replaced by D{t) — 
(1 — e~*/^i). This is our FD theorem: the fluctuations 
F are proportional to the dissipation D through a func- 
tion of the temperature. The fluctuations are absent at 
zero temperature, while they are maximal at inflnite tem- 
perature. Usually one considers the fluctuations of the 
one-qubit hamiltonian /i, in which case \Tv(haz)\ = 2E 
the splitting of the energy levels. 

We proceed now to discuss the link between dissipation 
and entanglement under two complementary standpoints 
(a third approach is given in Ref. |l^]). The last equality 
in (^ means that Uz{(j),9,a) is equivalent to Vz{(t>,9) 
up to local unitaries (LU); thus we can restrict to the 
Vz{(j), 9) for the study of entanglement. 

Dissipation and Entangling Power. The entangling 
power of a unitary operation U has been given different 
definitions ||l5|-[l7|. Here, we are interested in the cre- 
ation of entanglement during the thermalization process. 
In this context, the natural definition of the entangling 
power of a thermalizing machine V = Vz{4>, 9) is 



£[V] = ina.x£{Vp(E) , 



(20) 



with £(.) a measure of entanglement. As we said above, 
for two qubits there exist a measure of entanglement C, 
called concurrence, that is computable, basically by rank- 
ing the eigenvalues of a 4 x 4 matrix This may be a 
tedious task on the paper, but is a trivial one for a com- 
puter. Performing the optimization, we find for p > q 



C[Vz{<l),9)] ^C{VPi®£.Vl) ^psm2<p. 



(21) 



The maximal entanglement is thus produced when the in- 
put state is the excited state |1). The phases fluctuations. 
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measured by 0, don't show up: the entanghng power of 
Vz{4>,0) depends only on the dissipation, measured by 
(j). Moreover, if we want to introduce small fluctuations 
into the bath ^ (and a fortiori in the continuous-time 
limit), we must consider small values of 0. In this limit, 
the entangling power is increasing with dissipation. 

Equivalence under LU. We have noticed above that 
C/z (</!), 0, a) is LU-equivalent to Vz{(j),9). This is a par- 
ticular case of a general theorem stating that any 
unitary operation on two qubits is LU-equivalent to 



U, = 



where 



_ _ 4> 

case fix — — "2 



In our 
Since 



= y. 

e [0,f]; and fi, = -f 
<Jz] Vz ((/>, 9) , one can always choose 
9 G [— fi f ] within LU-equivalence. Now, for parame- 
ters in these ranges, the /i^ are uniquely determined p8| . 
Therefore Uz{4>, 9, a) is LU-equivalent to Vz{(f>' , 9') if and 
only if (/) = </)' G [0, ^] and 9' — 6'mod7r. In conclusion, 
two thermalizing machines are LU-equivalent if and only 
if they are dynamically equivalent. 

Irreversibility. We conclude on some general consider- 
ations about irreversibility. The thermalization process 
that we described is certainly reversible, since it involves 
only unitary operations. Thermalization appears as a 
consequence of entanglement, not of measurements. The 
environment does not play the role of measuring appa- 
ratus, but of "waste basket for information". In fact, 
the information encoded in the initial state of the sys- 
tem p is not lost after the thermalization process, but 
is encoded in a different way, being spread between the 
system and the bath. The initial state of the system can 
be reconstructed only if one knows which qubits of the 
bath have interacted with the system, and in which or- 
der. Without this knowledge, any attempt of reconstruc- 
tion of the initial state will fail [^. Thus, irreversibility 
arises here as the interplay of two information processes: 
(i) the quantum information on the initial state of the 
system is spread between the system and the bath, still 
in a reversible way; (ii) the classical information about 
the order of the collisions is lost, leading to the practical 
impossibility of running the process backwards. As an 
application, one can define "safes" for quantum informa- 
tion that can be "opened" with classical keys 

In summary: we have discussed the family of the ther- 
malizing machines for two qubits. These unitary opera- 
tions can be decomposed into two processes, dissipation 
and decoherence. Dissipation is linked to fluctuations — 
and this is expected, although our FD theorem is derived 
in a different framework than usual — and to the entan- 
gling power of the machine in the process. Both dissi- 
pation and decoherence are related to equivalence under 
local unitaries. 
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